Search for: All records

A three-dimensional, regional simulation is used to investigate ionospheric plasma density irregularities associated with Equatorial Spread F. This simulation is first driven with background electric fields derived from ISR observations. Next, the simulation is driven with electric fields taken from the WAM-IPE global model. The discrepancies between the two electric fields, particularly in the evening prereversal enhancement, produce disagreeing simulation results. The WAM-IPE electric fields are then studied through a simple sensitivity analysis of a field-line integrated electrodynamics model similar to the one used in WAM-IPE. This analysis suggests there is no simple tuning of ion composition or neutral winds that accurately reproduce ISR-observed electric fields on a day-to-day basis. Additionally, the persistency of the prereversal enhancement structure over time is studied and compared to measurements from the ICON satellite. These results suggest that WAM-IPE electric fields generally have a shorter and more variable correlation time than those measured by ICON. 

Abstract Continuous wave signals from a network of high frequency (HF) beacons in Peru and other instruments are used to reconstruct the regional ionospheric electron number density in the volume surrounding the network. The continuous wave (CW) HF signals employ binary phase codes with pseudorandom noise (PRN) encoding, and the observables include propagation time or pseudorange, Doppler shift or beat carrier phase, and amplitude. A forward model based on geometric optics in an inhomogeneous, anisotropic, lossy plasma is used to relate plasma number density to the observables. Plasma number density is parametrized in terms of a modified Chapman profile in the vertical and biquintic B‐splines in the horizontal. Sensitivity analysis is required both to model the ray amplitudes and to solve the two‐point boundary problem for each ray. Sensitivity analysis is performed here using reverse‐mode automatic differentiation. In particular, we use an LLVM compiler (Clang), the corresponding OpenMP library, and the Enzyme Automatic Differentiation Framework plugin to compute the sensitivity (gradients) of ray endpoints with respect to their initial bearings. The resulting algorithm exhibits no performance penalty compared to variational sensitivity analysis and is far simpler to implement. 

Abstract. It is generally accepted that modeling Farley–Buneman instabilities requires resolving ion Landau damping to reproduce experimentally observed features. Particle-in-cell (PIC) simulations have been able to reproduce most of these but at a computational cost that severely affects their scalability. This limitation hinders the study of non-local phenomena that require three dimensions or coupling with larger-scale processes. We argue that a form of the five-moment fluid system can recreate several qualitative aspects of Farley–Buneman dynamics such as density and phase speed saturation, wave turning, and heating. Unexpectedly, these features are still reproduced even without using artificial viscosity to capture Landau damping. Comparing the proposed fluid models and a PIC implementation shows good qualitative agreement. 

Abstract Observations of 30‐MHz coherent backscatter from sporadic‐Eionization layers were obtained with a VHF imaging radar located in Ithaca, New York. The volume probed by the radar lies at relatively high magnetic latitudes, on the northern edge of the mid‐latitude region and underneath the ionospheric trough. Banded, quasi‐periodic (QP) echoes observed from Ithaca are similar to those found in lower midlatitude regions. The Doppler shifts observed are smaller and, so far, do not appear to reach the threshold for Farley‐Buneman instability. However, many of the echoes exhibit fine‐scale structure, with secondary bands or braids oriented obliquely to the primary bands. Secondary bands have been seen only rarely at lower middle latitudes. In previous observations, the QP scattering has been linked to unstable neutral wind shears. Neutral wind shear commonly found in the lower thermosphere could play a key role in the formation of these irregularities and explain some morphological features of the resulting plasma density irregularities and the radar echoes. We consider whether neutral instability and turbulence in the lower thermosphere is the likely cause for some of the structuring in the sporadic‐Elayers. Results of 3D numerical simulations of atmospheric dynamics in the mesosphere to lower thermosphere support the proposition. In particular, we focus on Ekman‐type instabilities that, like the more common Kelvin‐Helmholtz instabilities, are inflection point instabilities, although specifically associated with turning shears, and result in convective rolls aligned close to the mean wind direction, with smaller‐scale secondary waves aligned normal to the primary structures.